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Abstract 

The photon is explored as an object to test the applications of QCD to the pertur- 
batively calculable collinear parton distributions. We investigate analytic properties 
of DVCS amplitudes and related sum rules of generalized parton distributions of the 
photon using as an input their earlier calculations in the leading order. The relation of 
these GPDs to the quintessential functions in the framework of the dual parametriza- 
tion approach is also found. 

1 Introduction 

In the concept of QCD factorization the amplitude of a hard process is expressed in 
terms of convolution of a hard perturbative amplitude and nonperturbative function 
characterizing nonperturbative strong interactions. In particular Generalized Parton 
Distributions (GPDs) and Generalized Distribution Amplitudes (GDAs) [ - ] specify- 
ing low momentum scale strong interactions are used in the exclusive processes. In this 
Letter, we explore the application of the QCD tools to the photon. The unique feature 
of the photon is that in the leading order its partonic content can be calculated in per- 
turbative QED [ , ] . Investigation of the partonic content of the photon (including the 
gluonic corrections) started with the seminal paper by Witten [ ]. More recently the 
deeply virtual Compton scattering (DVCS) amplitude with the photon target (Fig.l) 
was calculated, and then GPD and GDA functions characterizing quark content of 
the photon were extracted. DVCS on a photon is therefore a kind of toy model for 
QCD which has been recently generalized for the case of impact parameter dependent 
GPDs [7,8]. This provides an incredible opportunity to test our current mathematical 
tools of GPDs studies. Exploring the connection [9] of usual GPDs, where skewness 
1^1 < 1, with the Generalized Distribution Amplitudes (GDAa), where |^| > 1, we 
checked holographic sum rules [ , ] and calculated related D-terms [ ]. In this 
way GPDs became known in the full region of ^, and this allowed us to apply the 
inverse Radon transform (derived for GPDs in [ ]) to obtain double distribution (DD) 
functions. Later on, we investigated subtleties of applying the recently suggested pre- 
scription [ ] to photon DDs. We also turned to dual parametrization [ ], where 
the partonic structure is represented as an infinite series of t-channel exchanges, and 
derived the quintessence functions. 
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Figure 1: 7*7 77 scattering 

2 Photon GPDs in physical and nonphysical re- 
gions 

DVCS amplitude can be represented as the tensorial decomposition [15] 

T'^-/^(Ar = 0) = Igl^^gfA^+l (^r^T^ + 9T9f " 9^' 9^) A2+\ {9't9t^ - 9^9^) As 

(1) 

Amplitudes are presented as convolutions: 

1 1 
Alio = / dxCv{x,OHi{x,^,0), AsiO = / dxCA{x,OH3{x,^,0) (2) 



with 



(3) 



Dropping the factors ^^^^ In ^ (as we are interested in x,^ dependencies only), 
photon GPDs corresponding to respective amplitudes in (1) are [ ]: 



Hf{x,C,0)^9{x-0 



x' + {l-xf-e 



+ 



- x)9{x + - oi-^ - 



1-e 

x2 + (l + x)2-^2 



(4) 



Hlix, 0) = e{x - 



GDAs of the photon are [ ] : 



x^-{i+ x)^ - e 



(5) 



^liz', C, 0) = e{z' - C f'^^''^ ^'^ + 9{z' - O^^^^^r^ + 



(6) 



mz', c', 0) = 0{z' - C')^ - 0{z' - C)^ - 0{C' + eiC - z'f-^ (7) 
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For convenience we used a more symmetric way [ , ] to define the coordinates of 
the photon distribution amphtudes by the difference of their momenta: 



l-z 

/_ i-c 

2 ' 



(8) 



where z^ - are variables used in [5]. 

Photon GPDs obtained from the DVCS process are defined in the region with 
1^1 < 1 and \x\ < 1. They can be extended [ ] to the unphysical region of |^| > 1, 
expressing GPD via GDA ^(z^) in their physical region — 1 < z < 1, — 1 < < 1. 

In adopted normalization [2,4,5] the relation between GPD and GDA takes the 
form: 

H{x,^)^^-sgn{my^) (9) 
Applying (9) to (6) and (7) we get GPD Hf{x,^) for |^| > 1, 



(x-g)(2xg+g2-l 



-1 < X < 1 
X > 1 



(10) 



and analogously for H^{x,^), 



-1) 



2(x+g) 



X < 



-l< X <l 
X > 1 

X < -1 



(11) 



As a result the functions are known in the full definition region (|^|<1,— l<x<l 
and 1^1 > 1, — ^ < X < ^). Their 3D figures are illustrated in Figs. (2) and (3). Their ^ 




Figure 2: Hf{x,0 



Figure 3: if|(x,0 



slices are illustrated in Figs. (4) and (5) respectively. As one can see the x-derivatives 
of i7i(x,^,0) and Hs{x,^^0) have discontinuities at the points —1, — ^, ^, 1. 
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Figure 4: /^^^(x,^ for C = 0.2 (dashed), Figure 5: for ^ = 0.2 (dashed), 

l.l(sohd), 1.6 (dotted) l.l(sohd), 1.6 (dotted) 

3 Holographic sum rules 

Here we investigate dispersion relations connecting real and imaginary parts of the 
DVCS amplitudes [ , ] in apphcation to photon. The GPD contribution to hard 
exclusive amplitude (particularly for DVCS like in (2)) in the leading order is defined 
(dropping the — 2e^) through 



^i,3(C,^)= / dxHi^s{x,^,t)[^-- 



ie X — ^ + ie 



(12) 



or simply 



A 



1,3 



ImAi^six, t) = -7ri?i,3(x, x, t). (14) 

Note that the imaginary parts of the ampHtudes Ai and A^, in the region a: > 
differ only by the sign because 



Hi{x,x) — —H3{x,x) 



1 — X 

1 + x 



(15) 



The holographic sum rule for GPDs ['0,11] guarantees that the relevant informa- 
tion of the function up to subtraction is contained in the 1-dimensional section: 



1 
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H{x, ^) — H{x, x) 
X — ^ 



dx = A 



/ 



da 



a — 1 



(16) 



-1+ 



where P J means principal value of integral, while A does not depend on ^, as a result, 



„ P f ImA(x) , 

7T J X — q 



(17) 



Here and further we omit indices 1,3 for amplitudes, GPDs and DDs when the equa- 
tions are valid for both Hi and H^. This equation tells us that the DVCS amplitude 
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is defined through its imaginary part up to the subtraction constant. Let us stress 
once more, the information about DVCS in the leading order is contained on the hne 
X = manifesting its holographic property. 

To define G one should consider GPDs connection to Double Distributions through 
integral [ ] 



l-\a\ 



H{z, i)^ j da J dl3F{l3, a)S{z - (3 - a^) 



(18) 



-i+|a| 



which can be interpreted [ ] as the 2-dimensional Radon transform (RT) [ ] (being 
the basis of tomography methods, for a recent review see e.g. [ ]), an integral of 
F(/3, a) over the line with the slope —1/^ and crossing the /3 axis at the point x. See 
Fig. 6. The application inverse RT for GPD leads [ ] to 





a 




H(x,^) for ^>1 




\ \ ^^^^^^ 
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Figure 6: H{x,^) as integral of F{/3,a) 
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= /§/ daHiz + (3 + atO-H{P + a^,0) (19) 



-oo — oo 



There is an addition to the DD, the so called D-term [9, 12], introduced in order 
to preserve the polynomiality condition which resides \x\ < |^|, 

1 l-\a\ 

H{z,0 = Jda J d/3{F{/3,a) + ^6{/3)D{a))S{z- P-a^. (20) 
More generally [9] 

H{z,i)^ j da j dl3{F{l3,a)+^G{/3,a))5{z- (21) 



-1 -l+|a| 



l-|a| 



-l+|a| 



for the two-DD representation, where F and G are defined ambiguously and can be 
reduced to the form [ 1 



F{/3,a) = /3f{/3,a) 
G{l3,a) = af{/3,a), 



where / is defined in the single-DD representation: 

1 l-|a| 
— I atx 

-l+|a| 



1 ±-|u:| 

J da J d/3f{/3,a)S{ 



z — P — a^,) 



(22) 
(23) 



(24) 



Let us find the D-terms of photon GPDs(i7i and H3) exploring [ ] that 

D{a) = $(a,0), 

as a result one gets 

Di{a) = (|q;| - 1)(2|q;| + l)s5n(a) 
Dsia) = 



(25) 



(26) 



The plot of Di{a) in Figs. (10). It is interesting the sign of the D-term is formally 
in accordance with stability criteria for nucleons in vacuum [ ] , in nuclear matter [20] 
and Q-balls [21]. 

Let us now check the holographic sum rules for photon GPDs (16). For |^| < 1 
integration interval divides into 3 intervals 



/ 



H{x,0 - H{x,x] 



dx 



-1 -1 e 

which for Hi and i^a respectively gives: 

^-21og(e + l)-l + log(4) 



X — ^ 



dx. 



(27) 
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-H3(x 


x) 
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x) 
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-H3(x 


x) 


X 







c - 1 

g(g(log(4) - 2) + 4) - 4 log(g + 1) - 2 + log(4) 
^2-1 

3g-2 1og(g + l)-l+log(4) 
^ + 1 

2 (e log(e + 1) - ^(1 + log(2)) + log (i±i) + l) 



(28) 



4^1og(^) 



^2-1 

2 1og(g + l) + 2-log(4) 
C + 1 



Contributions from these intervals exhibit dependence on ^ as illustrated in 
Figs. 7,8, but their sum for ^ does not in accordance with holographic sum 
rule. 



Hi{x^ ^) — i/i(x, x) 



x 



dx = 21n2 for |{| < 1 



(29) 
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X — ^ 



dx = for 1^1 < 1 



(30) 



The calculated values are in accordance with (16) with r.h.s. (26). It is 
instructive to explore these integrals in detail. Expression in l.h.s. of (16) has 
a nontrivial limit <^ ^ 0. It happens that (cf. [ ]) this limit cannot be derived 
through substitution = in (16), because the central region [— C^C] becomes 
strictly zero and so does the central integral while its limit (28) does not: 



limP / ^i(^^O-^i(^.^) ^^^2-21n2^0, 



(31) 



what leads to discontinuity at ^ = 0. 





Figure 7: Contributions from integrals of Hi: J - 

-1 

^ 1 
red solid; J - green dashed; J - blue dotted 



Figure 8: Contributions from integrals of H^: J 

-1 

^ 1 
red solid; J - green dashed; / - blue dotted 



In full analogy convolution integrals for Hi and satisfy holographic sum 
rule in the region \^\ > 1 as well [19], 



i7i(x, ^) — Hi{x^ x) 
X — ^ 



dx = 21n2 for |^| > 1 



(32) 



Hsjx,^) - Hs{x,x) 
x-^ 



dx = for |<^| > 1 



(33) 



As we can see, subtraction constants do not depend on [ ], for any small 
but finite ^. 
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4 Recovering photon double distributions 

Different pairs F{/3^a) and used in the two-DD representation (21) are 

connected through gauge transformations [ ] , 

F(/3,a)^F(/3,a) + ^^^ 

G'(/3,a)^G'(/3,a)-^^^ (34) 

Let us find the one transforming the D - term type DDs to a) = /3f{/3^ a) 
and G{/3^a) = af{/3^a) appropriate to the single-DD representation [ ], 

F{/3,a) = l3f{l3,a) ^(/3^a) Fn{/3,a) 

G{P,a) = af{l3,a) ^ D{a) 

Fd corresponds to the (F^G) pair, where = 6{/3)D{a). Hd{x^C) is 

the part of GPD corresponding to the D-term. 

Subtracting Hd{x^^) = sgn{^)D{^) from the H[x^^) with the D-terms 
which we evaluated exphcitly (26), we proceeded with calculation of Fii)(/5,Qf), 
^3d(/5,o^) and /i(/5,a). Here we needed the inversion only for the single first 
term in the r.h.s. of (20) which is much simpler. Firstly we made inverse Radon 
transform numerically, then conjectured analytic ansatz for DDs, and finally 
checked it analytically. 

Fiz,(/3, a) = [2(1 - - |a|) - 1 + 5{a)]sgn{^), (35) 

which is shown in Fig. (9). 

Applying inverse Radon transform for H^[x^^) we obtain 

F:,D{^,a)=8{a)-l. (36) 

Hence we derive /i(/3, a) using the inverse RT for ^^^^'^"^ where one can identify 
the term 5(^)^ 

/i(/3, «) = ^ - 1 + 25{m - W\) + (37) 

IpI « 

This leads to 

Fi (/3, a) = 5{a)sgn{/3) - /3 (38) 



Gi(/5, a) = -a + 25(/5)(l - \a\)a + 5{/3)Di{a) (39) 
Now it is simple to derive the gauge transformation function (34) 

Xi (/?,«) = asgn{/3){l - \/3\ - \a\) + C (40) 
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Figure 9: Regular part of double distribution 




Since Xi is determined up to constant, we can set constant C to zero to make Xi 
vanish at the border of the support rhombus + |a| = 1, after it Xi is shown on 
the fig (11). Note that Xi can also be calculated using only the Gi function [ ] 

^ i-|a| 
= < 0) J dtGi{t,a)-e{f3>0) J dtGi{t,a) (41) 

Although D^{a) = 0, one may still perform the gauge transformation and write 
down the single-DD representation generating fictitious (genuine but not neces- 
sary because it can be completely eliminated by gauge transform) DD G3, 

n f(R ^ dxsila) ^^^^ 
- afsiP, a) = — 



aFsnip, a) = a — h /3 — (43) 
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X3(/3,«) 



-a)F3]j{t/3, ta)dt 



(44) 



X3(^,«) = «(1 



+ i«r 



shown in Fig. (11). Consequently, 



/3(/5,«) = 



8{a) sgn{/3) 



/3 (|/9| + |«|)2 



(45) 



(46) 





Figure 11: Gauge transformation functions respectively Xi(/9,tt) and X3(/9,tt) 



Now let us consider the plus- prescription notation [ ] for single-DD. Plus- 
prescription is a standard method for regularization of singularities for parton 
densities. For GPD it can be used to select part not containing D-term. The 
aim of this construction is to make /9-integrals of /(/5, a) finite. Note that plus- 
prescripted function is a distribution which is defined by equation: 



i-H 



[f{P,a)]^{K{P,a))d^= J f{P,a){K{^,a)-K{0,a))d^ (47) 



-l+|a| 



-l+|a| 



More formally, the general definition of DD " plus" term (presented by equation 
(75) in [1- ]) is: 



[/(/9,«)]+ = /(/5,«)-W J f{%a)d^. 

-l+\a\ 



(48) 



Surely the integral in the r.h.s. diverges (its divergence is in fact the reason 
to use plus-prescription). We will show that the application of this definition 
leading to equation (74) in [ 



fiP,a) = [fiP,a)],+SiP) 



D{a) 



a 



(49) 
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is incorrect in case of photon single-DD due to presence 5 {a). 
For our known a) general definition leads to 



LI/51 



l + 25(/5)(l- |«|) 



(50) 



At the same time from the single-DD approach it follows: 



(51) 



Verifying it for photon /i(/5, a) we see that it is correct. However, dividing it by 
a we will come to 

l-\a\ 



a 



j fii/3,a)dp 



(52) 



-l+|a| 



because of the presence of 5(a) in fi{/3, a) (37). Indeed, its integration gives 

l-|a| l-|a| 



Di{a) 



a 



+ 6{a) J ^ (53) 



-i+|a| 



-i+|a| 



where infinite discrepancy arises from the term containing S{a)^ which does not 
contribute in (51). 



5 Quintessence function 



An alternative way to introduce GPDs instead of Double Distributions is the 
dual parametrization which represents parton distributions as an infinite series 
of t-channel exchanges [22]. In the dual parametrization the amplitude is defined 
through the integrals [14]: 

1 

f dr 1 
W(e,t) = / -Nix,t)[^^^=] (54) 



i-Vi-g2 



ReAi^,t) = I —Nix,t)[-^ ^ + 



1 



/dr 1 2 

- (x, t)[ - ^==] + 2D{t) (55) 



'1 + f + x' 

which can be uniquely inverted in the following way 

1 
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Its physics content is more clear for Mellin moments of A^(x, t) which are directly 
related to the contributions of quark-antiquark states with definite angular mo- 
mentum in the t-channel [ 1: 



i 

/ 



dxx'^~'N{x,t) = - 



dz 



1-z ' 



(57) 



where is the distribution amplitude of exchange with the angular mo- 

mentum J as illustrated in Fig. 12. Substituting 7tH{^^ ^) we obtain quintessence 



Figure 12: GPD decomposition on t-channel exchanges with momentum J 



for the photon 



Ni{x) = -Nr,{x) 



{x — l){x'^ + 2a;log(a;) — 1) 
(x + l)2 



(58) 



illustrated in Fig. (13). As we see, A^i(x) = — A^3(x) because likewise do the 
imaginary parts of the respective amplitudes, see (see 15) 




Figure 13: Quintessence A^i(x) 



Equation (54) (after substitution ^ = |(1 + |) which is connected with 
Joukowskii conformal map) is exactly the Abel transform [ ] . This means that 
it can be interpreted as a slice of 2-dimensional axially symmetric function 



m{uj) 



M{uj) 



(59) 
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Figure 14: m{(jj) 




Figure 15: M{uj) 



For the photon 



ImA{^,t) = [ —M{uj,t)- 
J ^ 



m{u) = 



(60) 



(61) 



UJ{UJ + 1)3/2 

illustrated in Fig. 14, and M{u) = um{u) in Fig. 15 

6 Conclusion 

Dispersion relations for hard exclusive amplitudes are highly useful tools for 
QCD. In this sense, it seems extremely important to verify them. However, it is 
a rather difficult to test them for hadrons. That is why the photon is the unique 
object to analytically check the dispersion relations and related QCD tools. 
Also, the expressions of photon GPDs can give a hint for constructing hadron 
GPDs. At the same time the possibility of a direct experimental examination 
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of photon GPDs should not be ruled out. One of the processes which seems to 
be a possible application is the photon splitting in the nucleus electromagnetic 
field (see e.g. [ ]) when the outgoing photons are collinear to the initial. 

In this work we extended the original photon GPDs to the unphysical re- 
gion of |(^| > 1 using their GDAs. Having completely defined GPDs we verified 
the holographic sum rule (16) for the DVCS amplitudes in the leading order 
(17). We calculated the D-terms (26) and subtraction constants (29), (30). Also, 
we derived photon DDs for the single-DD (37), (46) and two-DD (35), (36) ap- 
proaches with the help of the inverse Radon transform. Afterwards, we verified 
equations for the plus-prescription and noted that the various definitions of the 
D-term via the single-DD function are not equivalent. Later, using the dual 
parametrization approach we calculated quintessence functions (58) for photon 
GPDs. As an outlook let us mention the possibility to generalize the suggested 
approach to impact parameter dependent Photon GPDs [ , ]. 
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